Entanglement spectrum in one-dimensional systems 
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We derive the distribution of eigenvalues of the reduced density matrix of a block of length I in a 
one-dimensional system in the scaling regime. The resulting "entanglement spectrum" is described 
by a universal scaling function depending only on the central charge of the underlying conformal 
field theory. This prediction is checked against exact results for the XX chain. We also show how 
the entanglement gap closes when £ is large. 

PACS numbers: 



The interest in quantifying the entanglement in ex- 
tended quantum systems has been growing in recent 
times at an impressive rate, mainly because of its ability 
to detect the scaling behavior in proximity of quantum 
critical points (see e.g. Refs. [1] as reviews). Among 
the various measures, the so-called entanglement entropy 
has by far been the most studied. By partitioning an 
extended quantum system into two blocks, the entangle- 
ment entropy is defined as the von Neumann entropy of 
the reduced density matrix pa of one of the two blocks. 
The success of this quantity can be understood because 
it is a single number able to capture the main features 
of the scaling behavior. In fact, in one-dimensional (ID) 
critical ground-states, when the block A is a segment of 
length t in an infinite system, the entanglement entropy 
diverges with the logarithm of the block size as [2-4] 
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where c is the central charge of the associated confor- 
mal field theory (CFT) and c[ a non-universal constant. 
Away from the critical point, Sa saturates to a constant 
value [3] proportional to the logarithm of the correlation 
length [4], 

However, the reduced density matrix (at least in prin- 
ciple) contains more information than the entanglement 
entropy. This information should be encoded in the full 
spectrum of the reduced density matrix, which we shortly 
call "entanglement spectrum" , following Ref. [5] . In this 
letter we calculate the entanglement spectrum for ID sys- 
tems in the scaling regime, i.e. at or close to a quantum 
critical point. The study of the distribution of eigenval- 
ues provides clearly a deeper theoretical understanding 
of entanglement and correlations in extended systems, 
but not only. Indeed, all numerical algorithms based on 
the so-called matrix product states (among which den- 
sity matrix renormalization group [6] is the best known) 
give as first output a truncated spectrum of the reduced 
density matrix and from this all the other quantities are 
derived. Consequently, the knowledge of some scaling 
properties of this distribution provides an optimal check 
for the convergence and the accuracy of the numerics 
and could be used for putting accurate bounds on the 
efficiency of these methods as already done from other 
quantities in Refs. [7]. 



The scaling behavior of the entanglement spectrum can 
be related to the properties of the moments of the re- 
duced density matrix: R a = Tr//| = X^A", where A-j 
are the eigenvalues of pa- In the scaling regime, R a can 
be written as 
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where c a is a non-universal constant and L e ff is the rel- 
evant length in the considered regime. For example, L e ff 
equals the length block £/a if A is part of an infinite gap- 
less system [2, 4] , L g — ^ sin ^ if A is in a finite gapless 
system of length L, and L e g = £/<z [4] if the system is 
gapped [valid when the correlation length £ (the inverse 
mass gap) is large, but smaller than all the other lengths 
like I, L] . Everywhere a stands for the scale setting the 
microscopic length, e.g. the lattice spacing. The same 
a dependence is found in the case of open systems (the 
exponent of L e g in Eq. (2) halves and so does Sa [4]) 
and when A consists of more disjoint intervals [4] and 
also in some non-equilibrium situations [8] . For practical 
reasons, it is convenient to write 



Ra 



', with 6= - lni cff > 0, (3) 
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where by simple inspection b is related to the maximum 
eigenvalue: b = — In A max , as well known [9]. 

In order to characterize the entanglement spectrum, 
we define the distribution of eigenvalues P(X) = S(X— 
Xi), which is normalized to m , where m is the dimen- 
sion of the local Hilbert space, e.g. m — 2 for spin 1/2 
systems. We determine this distribution for a ID system 
in the scaling regime where Eq. (2) applies. The only as- 
sumption in what follows is that the cv-dependence of the 
non-universal part of c a can be ignored, i.e. that we can 
write c a — a c ( Q_1 / Q )/ 6 f a , with f a a constant function 
(i.e. the effect of c a is to replace the lengths appearing 
in R a with the dimensionless quantity L e ff)- A priori 
this can appear as a very crude approximation, but it 
is not the case. For the gapless XX chain, c a is known 
analytically [10], and it is easy to show that f a varies 
less then 1% as soon as a > 2. Physically, this is equiv- 
alent to stating that the main contribution of c a is to 
set the microscopic length a to be used in the continuum 
description in terms of the lattice one. We will check a 
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posteriori for the XX chain that the result obtained ig- 
noring c a describes accurately the actual entanglement 
spectrum when L c ff is large. A very precise numerical 
determination of c a for the XXZ chain with — 1 < A < 1 
[11] shows that the same scenario is true on the full criti- 
cal line of the model. In few other cases c a is also known 
[12-14] and similar properties can be found. 

Thus, ignoring temporarily c a , we will compute the en- 
tanglement spectrum, using the simple observation that 
AP(A) = limIm/(A-ze), with 



If>„z-« = I / 



dX 
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Here, /(A — ie) has an imaginary part when e — > only 
on the support of P(A), due to the pole in the r.h.s. of 
(4). The calculation of f(z) is straightforward: 
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where Life(y) is the poly logarithm function, which is ana- 
lytic on the complex plane, with a cut on the real axis for 
y > 1 (that once again is just b = — lnA max ). The dis- 
continuity along the cut is given by lim ImLife(y + ie) = 

7r(ln7/) fc_1 /r(fc) for k > 1. The sum can be explictly 
done and we end up with 
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where the delta peak comes from the contribution of 
k = in (5) and Ik{x) stands for the modified Bessel 
function of the first kind. In this derivation, only R a 
with positive integers a enter. For these values (as we al- 
ready discussed) c a in general does not vary significantly. 
This gives an argument that explains why Eq. (6) works 
well for large enough L c g in the following calculations for 
the XX model and why the same is expected in general. 

P(e _t ) can also be obtained by considering the inverse 
Laplace transform of R a in the variable a. Using stan- 
dard results Eq. (6) is easily recovered. We preferred to 
give the previous derivation because it highlights the role 
played by positive integer a, where c a can be ignored. 
Conversely the inverse Laplace transform requires an in- 
tegral on the complex plane, over a contour where there 
is less control on the values of a that are contributing 
relevantly. We noticed this because numerical inverse 
Laplace transform could be done in principle in more 
complicated cases (like those in Refs. [10-14]) where the 
previous analytic reasoning does not work. 

Let us discuss now the main properties of P(A): 
(i) The mean number of eigenvalues larger than a given 
A is 



n(A) 



dAP(A) = 7 (2V^n(A max /A)) . (7) 



Note that for A — > 0, n(A) diverges, as it should, because 
in the continuum the number of eigenvalues is infinite. In 
the lattice models, this can be regularized by the finite 
number of degrees of freedom (as e.g. in spin chains), 
but not always (e.g. bosons have always infinitely many 
eigenvalues) . 

(ii) The normalization ^2 Aj = 1 corresponds to 
/ AP(A) = 1, and follows directly from Eq. (6). 

(iii) The entanglement entropy is given by 



S=- 



Aln XP(X)d\ 
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reproducing the result that the single copy entanglement 
equals one-half of the entanglement entropy [9]. 

(iv) Majorization is a relation between two probability 
distributions A = {Aj} and /j, = {^} whose elements 
are ordered Ai > A 2 • • • > A at (and similarly for y): 
it is said that A majorizes \i if J^tLi A« > Mi f° r 

any M = 1,...,N and £ti A; = ££i Mi = 1. It 
has been argued, observed numerically and in some in- 
stances proven analytically, that with increasing L e ff the 
resulting distribution of eigenvalues is majorized by the 
ones at smaller scaling lengths [15, 16] (sometimes this is 
referred to as majorization along renormalization group 
flow) . From the previous result it is straightforward that 
majorization holds in the scaling regime when Eq. (6) 
for P(A) applies. In fact, we have 



M 



S (M)ee]>> 



A„ 



1 



T 

Jo 



l (M) 



~ y2,Ab h{y) 



(9) 

which, at fixed M, is a monotonous function of A max 
(that is a monotonous function of L e ff)- This proves ma- 
jorization in a very easy way. It is also simple to check 
majorization directly by making the numerical integral 
at fixed M and varying A max . 

We now compare the scaling function for the entangle- 
ment spectrum with the eigenvalues of lattice models, in 
order to show its predictivity and to highlight its limits. 
As a prototype of lattice models we consider the gapless 
XX chain defined by the Hamiltonian 



(10) 



where a^' v are the Pauli matrices at site i. The reduced 
density matrix of a block of t contiguous spins in an infi- 
nite chain (so, among the cases of before, this corresponds 
to L e ff = I j a) can be obtained by exploiting the mapping 
to free fermions [3, 17]: 



Pi oc exp 



l<i,j<£ 



(11) 



where c, are spinless fermion annihilation operators, and 
* = M[(l-C7)/a C, = Sln ^:f 2 ■ (12) 
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FIG. 1: Sum of the first M eigenvalues of the XX model 
up to M = 100: 1 - s(M) as function of M for I = 
10, 100, 1000, 10000 (black dots). The red line is the confor- 
mal field theory prediction Eq. (9), in which A max has been 
fixed to the maximum eigenvalue obtained numerically. 



Calculating the eigenvalues z/, of the matrix dj re- 
quires only the diagonalization of an £ x £ matrix. In 
terms of the Vi, the 2 £ eigenvalues of the reduced density 
matrix are the products 



(13) 



where A is a subset of {1, 2, . . . , £} and B the comple- 
ment. While it is possible to obtain all the Vi up to very 
large £, making all of the 2 e products to obtain the full 
spectrum requires far too much memory on a personal 
computer. Thus we calculate the full eigenvalue spec- 
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FIG. 2: Inverse function of A; (value of the i-th eigenvalue 
of the reduced density matrix) for £ = 10, 100, 1000 for the 
XX chain. The plot is shown in terms of the scaling variable 
2^bln(Ama X /A). The full line is the CFT prediction. 
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FIG. 3: The same as in Fig. 2 for all 10 < £ < 28. Even (odd) 
£ are shown on the left (right) panel. 



trum only up to £ = 28 (which corresponds to almost 300 
millions of eigenvalues) , and for larger £ we truncated the 
fermionic spectrum to 24 modes (i.e. we diagonalize the 
C matrix finding all the z^, but in the products we do 
not vary the £ — 24 eigenvalues that are closer to or 
1 and generate the smaller A,). This reproduces in the 
exact order the first few thousands eigenvalues, to which 
we limit our analysis in the truncated cases. 

We start this analysis from the check of the "sum- 
rule" given by Eq. (9). In Fig. 1 we report 1 mi- 
nus the sum of the first M (up to 100) eigenvalues for 
£ = 10, 100, 1000, 10000. It is evident that increasing £ 
the sum is well described by Eq. (9) for larger and larger 
values of M. The (negative) deviations from the CFT 
take place at a value of M that roughly scales like ln£ 
and they are obviously due to lattice effects, because the 
sum rule must be saturated by a finite sum up to 2^ that 
cutoffs the integral of the continuum limit. 

A stronger test of our predictions is provided by Eq. 
(7) that gives the total number of eigenvalues between A 
and A max . This formula provides the natural scaling vari- 
able x = ln(A m ax/A). Once A max is determined from 
the numerics, n{X) is just given by Iq(x), independently 
from any other detail of the model under consideration 
and, more surprisingly, it is also independent from the 
model itself. Instead of plotting the number of eigenval- 
ues larger than a given A, we prefer to show the inverse 
function of \ that is the value of the i-th eigenvalue. In 
addition to represent n(A), this also has the advantage 
of giving information about the degeneracy of all eigen- 
values, as the number of points at the same position on 
the horizontal axis. For £ = 10, 100, 1000 this function is 
shown in Fig. 2 and compared with n(A) in Eq. (7). It 
is evident that when the spectrum becomes almost con- 
tinuous, it is well described by the CFT prediction up 
to a given A m i n (€) (i.e. a given x max {l)). Increasing £ 
the number of eigenvalues described by Eq. (7) increases 
and at £ = 1000, approximately thousands eigenvalues 
fall on the CFT scaling curve. The negative deviations 
for eigenvalues smaller than A m i n (i?) are due to lattice ef- 
fects as in the previous case of the partial sumrule. The 
degeneracies of the eigenvalues are not reproduced by our 
approach. In Fig. 3 we plot the entanglement spectrum 
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for all £ from 10 to 28 in order to show both finite-size 
and parity effects. When £ is odd (a case rarely consid- 
ered in the literature) the approach to the asymptotic 
scaling is much slower. This can be understood because 
of the double degeneracy of all the eigenvalues (included 
the largest), which moves a large weight toward x = 
that, because of the sumrule, must be compensated by a 
smaller weight at large x (i.e. small eigenvalues). 

There is a final interesting feature that is clearly vis- 
ible in Figs. 2 and 3. When plotting in terms of 
x = 2- v /61n(A max /A), the first few eigenvalues do not 
change their positions with changing I. The logarithms 
of these discrete eigenvalues result to be equispaced, i.e. 

61n^-fc^ ^ e -i^ ; ( M ) 
A y 

where [i and v are indices referring to two non-degenerate 
consecutive eigenvalues and k is a constant that is differ- 
ent for even and odd I (and we used b = — In A max = 
c/6ln£ and c = 1). When A M = A max , this simple for- 
mula tells us how the gap between the first and the sec- 
ond eigenvalues closes when increasing £, which, following 
Ref. [5], we call "entanglement gap". This scaling of the 
discrete part of the spectrum can be related to an old 
result about the scaling of the eigenvalues of the corner 
transfer matrix [18]. In fact, because &ln(A M /A„) is al- 
most independent of £, for this part of the spectrum the 
result should be the same as that of a segment £ in a 
system of length 2£. In this case the reduced transfer 
matrix can be seen as the fourth power of a corner trans- 
fer matrix [4] (of angle ir/2), for which the analogous of 
Eq. (14) has been derived from CFT [18]. Using this 
correspondence the asymptotic behavior of P(A <C 1) for 
gapped systems has been already derived in Ref. [19]. 
However our result goes far beyond, explaining also the 
reason of such a simple result. 

In conclusions, the main result of this letter is the an- 



alytic derivation of the universal entanglement spectrum 
for ID models in the scaling regime, given by Eq. (6). 
This turns out to depend only on the central charge of 
the underlying conformal field theory. We found that 
Eq. (6) describes accurately the continuum part of the 
spectrum of the XX chain for large enough L a g, and is 
expected to work for any model at or close to a quantum 
critical point. This can be seen as a surprise because 
it means that the continuum part of the entanglement 
spectrum does not contain more information than the en- 
tanglement entropy that is just one possible average over 
the spectrum. The reason for this (maybe unexpected) 
result can be traced back to the fact that conformal in- 
variance in ID is so strong that it fixes completely the 
shape of the full spectrum, leaving only one parameter 
(the central charge) free. Such parameter can be fixed 
by of a single "measure" like the entanglement entropy, 
the largest eigenvalue, etc. 

Oppositely, the discrete part of the spectrum is only 
reproduced in average by Eq. (6). As a consequence, 
the location and the degeneracies of the low-lying eigen- 
values of the reduced density matrix can still be a tool 
for extracting further universal information on the model 
under investigation that could not be captured by the 
entanglement entropy. For example, it is known that 
the degeneracies of eigenvalues of the isotropic Heisen- 
berg antiferromagnetic chain are larger than the XXZ 
ones (because of the larger symmetry) , but they all have 
c = 1 and so the same continuum part of the spectrum. 
A careful study of these issues is needed, but it goes far 
beyond the goals of this letter. 
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